function [] = fig_sec9_5_persitenceProbabilityAdjustmentPrimal()

close all;

ue = load ('data/persistenceProbAdjustmentPrimal_ue.txt');
pe = load ('data/persistenceProbAdjustmentPrimal_pe.txt');
objFunc = load ('data/persistenceProbAdjustmentPrimal_objFunc.txt');
jom_ue = load ('data/persistenceProbAdjustmentPrimal_jom_ue.txt');
jom_pe = load ('data/persistenceProbAdjustmentPrimal_jom_pe.txt');
jom_objFunc = load ('data/persistenceProbAdjustmentPrimal_jom_objFunc.txt');

multiplot ('Time', 'Network utility' , objFunc , jom_objFunc , [-900 -750] , '-' , 'p_e' , pe , jom_pe , [0 0.1] , '-' , 'u_e' , ue , jom_ue , [0 0.04] , '-');

% multiplot ('Time', 'Network utility' , objFunc , jom_objFunc , [-900 -500] , '-' , 'p_e' , pe , jom_pe , [0 0.2] , '-' , 'u_e' , ue , jom_ue , [0 0.1] , '-');
% multiplot ('Time', 'Network utility' , objFunc , jom_objFunc , [-100 100] , '-' , 'p_e' , pe , jom_pe , [0 1] , '-' , 'u_e' , ue , jom_ue , [0 1] , '-');

jom_ue
jom_pe
jom_objFunc
grad_ue = ue(end,2:end)
grad_pe = pe(end,2:end)
grad_objFunc = objFunc (end)
min(pe(end,2:end))
